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Abstract 


Experimental evidence exists which suggests turbulent boundary layer relaminarization 
may play an important role in the inverse Reynolds number effect in high-lift systems. An 
experimental investigation of turbulent boundary layer relaminarization has been undertaken 
at the University of Notre Dame’s Hessert Center for Aerospace Research in cooperation 
with NASA Dryden Flight Research Center. A wind tunnel facility has been constructed 
at the Hessert Center and relaminarization achieved. Preliminary evidence suggests the 
current predictive tools available are inadequate at determining the onset of relaminarization. 
In addition, an in-flight relaminarization experiment for the NASA Dryden FTF-II has 
been designed to explore relaminarization at Mach and Reynolds numbers more typical 
of commercial high-lift systems. 
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1 Introduction 


In a 1996 report entitled “Key Topics for High-Lift Research: A Joint Wind Tunnel/Flight 
Test Approach” [1] supported by NASA Dryden, Thomas et al. proposed a joint experimental 
investigation of a to-be-determined high-lift building block flow utilizing the University of 
Notre Dame’s Hessert Center’s low speed wind tunnels and NASA Dryden ’s Flight Test 
Fixture (FTF-II). Following the recommendation outlined in the report, an internal study 
conducted by the authors was initiated to determine which high-lift building block flow 
would be investigated. Results from an extensive literature search and high-lift research 
conducted at the University of Notre Dame, led to the selection of turbulent boundary 
layer relaminarization as the most suitable high-lift building block flow for the joint wind 
tunnel/flight test investigation. 

1.1 Relaminarization &; High-Lift Systems 

In a high-lift system, there is evidence to show^ that relaminarization of the main element 
boundary layer may have a profound effect on the global aerodynamic behavior of the system. 
If one examines the effect of Reynolds number on Ci mhX of a high-lift system, experience 
dictates that for increasing Reynolds number Ci mn would increase. In fact, this is not 
always the case as can be seen in Figure 1, w r hich demonstrates an example of Ci mAX behavior 
for a high-lift system. As one can see, the low Reynolds number region is approximately 
linear, but as Reynolds number increases Cl raax decreases dramatically until at even larger 
Reynolds number it begins to increase again. Although it is not illustrated in Figure 1, 
at even greater Reynolds numbers flight test experiments have shown that Ci mmx can drop 
off again. This phenomena is referred to as the “Inverse Reynolds Number Effect” and its 
cause is unknown. If one w*as interested in.high Reynolds number behavior, as is often the 
case for flight aircraft, merely extrapolating from low Reynolds number data typical of wind 
tunnel tests would be insufficient. 

Flight test experiments on leading edge transition and relaminarization conducted by van 
Dam et al. [2] using the NASA Transport Systems Research Vehicle, a Boeing 737-100, has 
provided tantalizing evidence, but not proof, that relaminarization may be responsible for the 
inverse Reynolds number effect that can occur at high Reynolds numbers. Relaminarization 
in high-lift systems may occur as the turbulent boundary layer proceeds from the attachment 
location, under the leading edge, around the nose of the main element. As the turbulent 
boundary layer proceeds around the nose, it encounters a strong favorable pressure gradient. 
The accelerated boundary layer w r ill then merge with the leading edge slat w y ake over the main 
element forming a confluent boundary layer. This thick viscous layer will then encounter an 
adverse pressure gradient as it moves downstream over the main element ultimately leading to 
flow separation, which has an unfavorable effect on Cl max - However if a sufficiently strong 
favorable pressure gradient exists at the leading edge of the main element, the turbulent 
boundary layer may relaminarize, thus thinning the boundary layer and resulting in a delay 
of confluence. Delayed confluence will move the separation point farther aft on the main 
element, wdiich has a favorable effect on Thus a possible failure of the turbulent 

boundary layer to relaminarize at high Reynolds number may be responsible for the inverse 
Reynolds number effect. 
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Figure 1: Taken from Meredith, P.T., “Viscous Phenomena Affecting High-Lift Systems and 
Suggestions for Future CFD Development, 11 AGARD CP-515, pp. 19-1, 19-8, 1993. 


The criterion for onset of relaminarization is typically attributed to the relaminarization 
parameter K\ a dimensionless group which results from a non-dimensionalization of the x- 
momentum Navier-Stokes equation and is comprised of the pressure gradient term, the 
freestream velocity, U^, and the kinematic viscosity, v y K ) as defined by Launder [3], is 
typically written as, 


K = v dU a c 
VI dx 


or K can be rewritten in terms of the pressure coefficient C p and the Reynolds number, Re, 
as 


K = -I-L( J 

2Re'l-C,' dx 


Low Reynolds numbers experiments have shown relaminarization typically occurs for values 
of K > 3 x 10“ 6 . Upon inspection of the relaminarization parameter, it is evident that 
K is determined solely from external inviscid flow parameters and thus can be estimated 
using only surface pressure measurements. However one would expect that boundary layer 
parameters ought to influence relaminarization. 

The presence of the — term in Equation 2 should not be surprising if the mechanism 
for the inverse Reynolds number effect truly is relaminarization. In effect, at 1ow t Reynolds 
numbers K is large ( > 3 x 1G -6 ) and if K is the indicator for relaminarization, one would 
expect the leading edge turbulent boundary layer to be in a relaminarized state. Equation 
2 indicates that as Reynolds number increases, K decreases and one would expect that once 
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the Reynolds number became large enough for the relaminarization parameter to achieve a 
value less than the critical value of 3 x 10 -6 , relaminarization would cease. Focusing our 
attention on a high-lift system, at low Reynolds number, one would expect the main element 
turbulent boundary layer, as it moves from the attachment line around the leading edge, to 
relaminarize. As Reynolds number increases, a critical Reynolds number will be reached. 
Above this critical Reynolds number, K will drop below the critical value of 3 x 10“ 1 2 3 4 5 6 , which 
will lead to the cessation of relaminarization. Cessation of relaminarization will then lead to 
earlier confluence of the main element boundary layer and slat wake, followed by a upstream 
relocation of the separation point on the main element, which will result in a decrease in 


1,2 University of Notre Dame/NASA Dryden Research Approach 

Experiments conducted in academic wind tunnels at low Reynolds numbers have shown that 
for K > 3 x 10~ 6 relaminarization will occur. However at large Reynolds numbers, typical of 
those found in flight, the likelihood of relaminarization occurring is an open question. The 
uncertainty surrounding the importance of relaminarization in high-lift systems and the lack 
of fundamental knowledge of the mechanisms responsible for boundary layer relaminarization 
have motivated the current joint wind tunnel/flight test project between NASA Dryden 
Flight Research Center and the University of Notre Dame. This joint effort involves wind 
tunnel testing at the University of Notre Dame in support of future flight tests on the NASA 
Dryden FTF-II. 

Two goals have been established for wind tunnel testing at the University of Notre Dame: 

1. To better understand the flow physics of turbulent boundary layer relaminarization, 
and 

2. To develop diagnostic tools in support of FTF-II flight testing at Mach and Reynolds 
numbers typical of commercial aircraft. 

The remainder of the report will discuss how these goals were achieved and to what 
degree. More specifically, the following will be discussed: 

1. Defining relaminarization and how best to experimentally determine when relaminar- 
ization is present, 

2. Describing the relaminarization wind tunnel facility and its capabilities, 

3. Documenting the experimental relaminarization results, 

4. Defining a possible relaminarization design concept for the NASA Dryden FTF-II, and 

5. Recommendations for future experimental measurements to further our knowledge of 
relaminarization and to aid in the development of the NASA Dryden FTF-II relami- 

narization experiment. 
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2 Relaminarization 


2.1 Background 

The first evidence of turbulent boundary layer relaminarization, or “reverse transition”, was 
documented by Wilson and Pope 1954[4] while studying heat transfer coefficients on gas 
turbine blades. In the nearly five decades since Wilson and Pope’s findings, there have 
been several important experimental investigations focusing on turbulent boundary layer 
relaminarization. Most notably Launder 1964[3], Patel and Head 1968[5], Blac.kwelder and 
Kovasznay 1972, Narasimha and Sreenivasan 1973[6],and more recently Fernholz and War- 
nack 1998[7],[8] and Ichimiya, et al. 1998 [9]. Also, an excellent review article was published 
by Sreenivasan [10] in 1982. The experimental investigations conducted by these researchers 
have provided much information regarding the processes and mechanisms associated with 
turbulent boundary layer relaminarization. But as will be shown, there still exists many 
unanswered questions surrounding the prediction of relaminarization and its occurrence out- 
side a laboratory environment. 

2.2 The Mechanics of Relaminarization 

Relaminarization, or “reverse transition”, occurs when a turbulent boundary layer is exposed 
to a large enough favorable pressure gradient over a suitable duration causing the turbulent 
boundary layer to revert to a laminar-like state[3],[6],[8],[10],[ll],[l2]. Relaminarization is 
characterized by: a thinning of the boundary layer, a departure of the mean velocity profiles 
from both the “law of the wall” and “law of the wake” profiles [5], an initial decrease then 
rapid increase in the shape factor [6] ,[11], an initial increase in the heat transfer and skin 
friction coefficients followed by a substantial decrease[5],[6],[ll], a decrease in the relative 
turbulence intensity [11], [13], [14], a rapid decline of turbulent bursting in the wall-layer [15], 
a spread of turbulent intermittency from the outer-layer to the wall-layer [13], and a decay 
of the turbulent stresses in the near-wall region [8], [14]. 

For a boundary layer undergoing relaminarization, Launder[3] documented an energy 
shift toward the low wave number spectrum in the near- wall region and a similar, more 
pronounced, shift in the wave number spectrum of the outer-region. Indeed it has been 
observed that the turbulent spectra of a relaminarizing boundary layer deviates from the 
typical roll-off associated with spectra of a fully developed turbulent boundary layer. 
The total energy associated with the relaminarized boundary layer is in fact smaller than 
the total energy of a fully developed turbulent boundary layer in a zero pressure gradient 
environment. 

For a turbulent boundary layer undergoing relaminarization, Fernholz and Warnack[8] 
found that the integral length scales, defined as 

A* = A r u (3) 

j ry+Ay 

A V = o / Ru'u'dAy (4) 

* Jy-Ay 
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where A r is defined as 


A t = f R r dr (5) 

Jo 

show an increase in A X /S by a factor of four and an increase in A y /6 by a factor of two 
over a typical turbulent boundary layer in zero pressure gradient. An elongation of the 
large eddy structure in the streamwise direction due to the large favorable pressure gradient 
is also evident in a relaminarized turbulent boundary layer. These elongated large scale 
structures show no evidence of scaling with the boundary layer thickness, <5, as they do in 
zero and adverse pressure gradients. When emerging from the strong favorable pressure 
gradient region into a zero pressure gradient environment, the large scale structure contain 
definite history effects attributed to flow acceleration. 

2,3 A Model For Relaminarization 

In describing the sequence of events that would lead to relaminarization, Launder [3] com- 
mented that as a turbulent boundary layer is exposed to a large favorable pressure gradient 
the turbulent structure in the near-wall region begins to depart from its fully turbulent state. 
The turbulent boundary layer becomes unable to respond to the large flow acceleration and 
a rapid momentum increase results. This increase in momentum causes an increase in the 
viscous shear stress, which if the pressure gradient persists, grows larger than the turbu- 
lent stress due to the inability of the turbulent structure to adjust to the flow acceleration. 
Consequently, the turbulence dissipation exceeds production, which leads to a decay of the 
turbulent structure, resulting in relaminarization. 

Elaborating on Launder’s premise, Sreenivasan[10] incorporated a two-layer model com- 
prised of a viscous inner-layer and an inviscid outer-layer to explain the mechanism responsi- 
ble for relaminarization. Sreenivasan argued that when a turbulent boundary layer is exposed 
to a large favorable pressure gradient, the turbulent structure is distorted in the outer-layer 
by the rapid flow acceleration. During this rapid acceleration, a new viscous inner-layer 
develops in which the initial turbulence decays. Unlike the inner-layer, turbulence decay in 
the outer-layer is insignificant. The two layers interact weakly, in so far as they only pro- 
vide the appropriate boundary conditions. The major contributor to relaminarization in this 
view is the “domination of pressure forces over the slowly responding Reynolds stresses in 
the outer-layer, accompanied by the generation of a new laminar subboundary layer, which 
itself is maintained stable by the acceleration.” 


2.4 A Criterion for Relaminarization 

Relaminarization is not a catastrophic phenomenon and, as such, the onset of relaminariza- 
tion is invariably difficult to define. In his 1982 review paper, Sreenivasan [10] assembled a 
number of the typical parameters that have been used to define the onset of relaminarization; 
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Parameter 

Numerical Value 

K 

~ 3.0 x 10" e 

Ap 

~ -0.025 

At 

~ 0.025, n = — | 

Kc n f 

~ 0.025, n = - 1 

H 

min ( H ) 


Table 1: Common relaminarization parameters and their values at the onset of relaminar- 
ization 


they include K , A p ) A r } Kc n ^ and H where these quantities are defined as 


K = 
A p = 

At = 

yjn 

xc f = 

H = 


v du 

u 2 dx 
v dP 


u* dx 



v du n 

u 2 dx f 


S' 


0 


(6) 

(7) 

( 8 ) 
(9) 

( 10 ) 


Typical values of the above parameters used as indicators of relaminarization are listed in 
Table L 

Narayanan and Ramjee[ll] observed the onset of relaminarization to correspond with 
a decrease in the turbulence intensity. They also observed the disappearance of turbulent 
bursts near the wall occurring for similar values of the above relaminarization parameters. In 
addition to the previously mentioned criteria } Fernholz and Warnack[8] recently documented 
two higher moments, the skewness, S Tw ,, and the flatness, F Twh of the fluctuating skin friction 
coefficient reaching maximums at the initiation of relaminarization. 

Although the above parameters are useful in approximately determining the onset of 
relaminarization, they are often not adequate enough and offer little to no help in explain- 
ing why a turbulent boundary layer will relaminarize. In fact, SreenivasanflO] comments 
“...these criteria are often useful as rough indicators of some drastic changes in the acceler- 
ated turbulent boundary layer.” And “...that the convenience associated with the use of any 
of these criteria should not mask our concern about their basic inadequacy as markers of the 
‘onset’ of relaminarization.” Also Narasimha[6] writes “...it would be somewhat surprising 
if, over an appreciable Reynolds number range, the relevant parameter were to be one like 
K which takes no account of the shear flow at all.” 
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3 Experimental Work Completed 

3.1 Relaminarization Test Facility 

To explore the role of relaminarization in high-lift systems and to gain a better under- 



Figure 2: Hessert Center relaminarization wind tunnel test facility. 


standing of the mechanism(s) responsible for relaminarization, a wind tunnel test facility 
was constructed at the Hessert Center utilizing the 0.6 m x0.6m (2 ft X 2 ft)subsonic indraft 
wind tunnels. The design of the w r ind tunnel test section was motivated by the desire to 
expose a turbulent boundary layer to a region of favorable pressure gradient corresponding 
to a constant K value. Inspection of Equation 2 reveals that at a fixed Reynolds number, 
to maintain a constant value of K over a specified streamwise distance, it is necessary that. 



(ii) 


where C is a constant. Invoking conservation of mass and the definition of C p it can be 
shown that the required tunnel geometry to maintain a constant value of A" is a linear profile 
given by, 


h(x*) = h Q (1 — K Re L x *) (12) 

where h is the height of the wind tunnel section, h 0 is the initial height of the converging 
section, Re*, is the Reynolds number based on L, the length of the converging section, and 
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Figure 3: Side view of the current relaminarization wind tunnel facility at the "University of 
Notre Dame. 



Figure 4: Internal view of the converging section of the relaminarization wind tunnel facility 
at the University of Notre Dame. 
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x * is the streamwise distance non-dimensionalized by L. Equation 12 implies that for a 
fixed tunnel geometry, changing the freestream velocity results in changing the value of K. 
Thus the wind tunnel section was designed to provide a constant K value over a distance 
of 2 ft for a range of velocities that would correspond to values of K between approximately 
5x 10 -6 and 1 x 10 -6 ; values large enough to obtain relaminarization and values small enough 
where one would expect relaminarization not to occur. Figure 2 provides a schematic of the 
current relaminarization test facility, while Figure 3 is a side photograph of the wind tunnel 
test section and Figure 4 is a photograph showing an internal view of the converging section. 


Comparison of Experimental C p Distribution with the Theoretical Prediction 



Figure 5: Theoretical prediction of C p distribution overlayed on experimental C p curve. 
Demonstrates a region of approximately 0.43 m (17 in) where K is constant. 


Figure 5 contains a plot of the theoretical pressure distribution for a constant value of K 
compared with the experimental pressure distribution measured in the converging section of 
the wind tunnel. The theoretical pressure distribution is represented by the solid curve in 
Figure 5 and is described by 


a 


i)- (i+Cj 


Pi =3. 35 m 


)- 


1 

1 — K Re/, 


(13) 


which satisfies the relation in Equation 11. The data points in Figure 5 represent the 
experimentally measured coefficient of pressure in the converging section of the wind tunnel. 
The theoretical and experimental pressure distributions are in acceptable agreement over a 
range of approximately 0.43 m (17 in) , w T hich corresponds to 


Ax 0.43 m 

L 0.61 m 


0.71 


(14) 
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Figure 6: Range of attainable K values for corresponding wind tunnel freestream velocities. 
K = 3 x 1CT 6 corresponds to Re L ~ 220,000. 

or 71% of the converging section of the wind tunnel (where L is the length of the converging 
section). This demonstrates that the facility has the ability to maintain a region of constant 
K over approximately this streamwise distance (0.43 m (17 in)). This length is found to 
be remarkably consistent over the designed operating range of the wind tunnel test section. 
Figure 6 illustrates the attainable K values for the designed operating range of the wind 
tunnel facility. 

3.2 Experimental Results from the Relaminarization Test Facility 

3.2.1 Initial Conditions 

To document relaminarization and determine the mechanisms responsible for relaminariza- 
tion, it is necessary to have a fully developed turbulent boundary layer. Thus the first tests 
conducted in the relaminarization facility involved documenting the development and state 
of the initial boundary layer upstream of the favorable pressure gradient region. Maint aining 
a turbulent boundary layer throughout the operating velocity range was more difficult than 
first thought. Many changes to the boundary layer development region were explored to en- 
sure a fully developed turbulent boundary layer w T as present upstream of the relaminarization 
region. 

Using a 90° X-wire boundary layer probe, a full flow field survey was conducted upstream 
of the relaminarization region to document the boundary layer state. Because the flat plate 
boundary layer develops in a nominally zero pressure gradient environment upstream of the 
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relaminarization region, the fluctuating velocity profiles were compared to Klebanoff’s[16] 
1955 data for a fully developed turbulent boundary layer on a flat plate. Figure 7 is a plot 
of the turbulence intensity versus the normal distance from the wall nondimensionalized by 
the local boundary layer thickness. Figure 8 is plot of the fluctuating v velocity scaled by 
the freest ream velocity plotted versus the normal distance from the wall nondimensionalized 
by the local boundary layer thickness. Figure 9 is plot of the Reynolds stress plotted versus 
the normal distance from the wall nondimensionalized by the local boundary layer thickness. 
In Figures 7-9, the profiles showm are for 4 < £/<*, < 9 ™ at x = 3.35 m (132 in) (the location 
directly upstream of the converging section of the wind tunnel) and Klebanoff’s profiles are 
also plotted for comparison. 

The turbulence intensity profiles (Figure 7) shows good agreement with Klebanoff’s data 
for 6 < U 0 © < 9 S. The = 4,5 j deviate from the other profiles for 0 < y/6 < 0.5. 
Likewise, the Reynolds stress profiles in Figure 9 display good agreement with Klebanoff’s 
data. The fluctuating normal velocity component profiles show good agreement with one 
another, but are shifted off Klebanoff’s curve; this is believed to be a bias error caused by 
alignment of the X-wire probe. 


Turbulent Spectra 

y/6 = 0.25 



K {m* 1 ) 


Figure 10: Turbulent u component velocity spectra comparison for V ^ — 4 — and U x = 8 — 
at x = 3.35 m (132 in) . 

The discrepancy between the measured profiles, specifically U 0 0 = 4 — , led to an ex- 
amination of the turbulent spectra. The ^-component velocity spectra for = 4 ^ and 
Uqo = 8 — is plotted in Figure 10 and clearly the turbulent spectra exhibit a region in which 
the roll-off is consistent with the expected K~% for turbulent boundary layers. Thus based 
on the behavior of the mean profiles and an examination of the turbulent spectra, it was 
concluded that for all initial velocities in the range of 4 “ < Uoo ^ 9 “ there was indeed a 
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u — 

'-'do s 

6 [mm] 

9 [mm] 

6* [mm 

H 

Re fl 

C f 

II ” 

4 

60.0 

6.08 

8.91 

1.47 

1545 

0.00390 

0.1639 

5 

57.5 

5.33 

7.41 

1.39 

1770 

0.00400 

0.2168 

6 

60.0 

5.86 

7.85 

1.34 

2344 

0.00380 

0.2548 

7 

60.0 

5.41 

7.22 

1.34 

2556 

0.00375 

0.2990 

8 

60.0 

5.31 

7.08 

1.33 

2839 

0.00370 

0.3359 

9 

60.0 

5.12 

6.85 

1.34 

3042 

0.00365 

0.3711 


Table 2: Summary of boundary layer parameters for the initial turbulent boundary layer at 
x ~ 3.35 m (132 in) 


initially turbulent boundary layer present. Table 2 contains a summary of several boundary 
layer parameters for the initial turbulent boundary layer at x = 3.35 m (132 in) . 



3.2,2 Relaminarization Results 

While testing was conducted for six discrete K values, results will be discussed for only two 
cases: 


1. K = 4.10 x ltr 6 , U 00 = i f 

2. K = 2.05 x 10’ 6 , U 0 0 = 8f 

For Case 1, we would expect to observe relaminarization since K > 3.0 x 10 6 , however 
for Case 2 we would not expect to see relaminarization since if < 3.0 x 1CT 6 . Examining 

Figure 11, a plot of the boundary layer thickness evolution, it is apparent that the Case 

1 boundary layer thins at a greater rate than the Case 2 boundary layer. The Case 1 
boundary layer also is thinner than the Case 2 boundary layer in the converging section 
although the Case 2 boundary layer is moving at twice the initial velocity; this is consistent 
with relaminarization occurring in the Case 1 boundary layer. Focusing on the shape factor 
evolution, where 



(15) 


inspection of Figure 12 reveals similar trends between the Case 1 and Case 2 boundary layers. 
At the onset of relaminarization, we would expect to see a sharp decline in the shape factor, 
in fact reaching a minimum value, followed by a gradual increase. Case 1 demonstrates 
this behavior in the region of x = 3.94 m (155 in), but it is not apparent to the author 
how to differentiate between the behavior in Case 1 and Case 2 as relaminarization and no 
relaminarization present. Examination of the skin friction values in Figure 13 shows a large 
decrease in the Case 1 boundary layer in the region of x = 3.94 m (155 in) when compared to 
the Case 2 boundary layer. This behavior is consistent with the onset of relaminarization in 
the Case 2 boundary layer. Figure 14 illustrates the evolution of the maximum turbulence 
intensity; it is evident that the Case 1 boundary layer sees a reduction and maintenance in 
turbulence intensity levels when compared to the Case 2 boundary layer, which continues to 
grow with downstream distance. This behavior would suggest relaminarization occurring 
in the Case 1 boundary layer. Finally, examining the log-law profiles in Figure 15, there is 
a noticeable deviation from the expected log-law behavior in the x = 3.94 m (155 in) plot; 
typical of a boundary layer undergoing relaminarization as observed by Patel and Head[5]. 
However, this behavior is observed for both Case 1 and 2 when one would expect only to 
see a deviation in the Case 1 boundary layer, if deviation of the boundary layer profile from 
the typical log-law profile is a proper criterion for the onset of relaminarization. In fact, 
Fernholz and Warnack[7] found similar deviations in log-law behavior for boundary layers in 
favorable pressure gradients with K < 3.0 x 10“ 6 and no observable relaminarization. 

From the above data, it is obvious to see that determining the onset of relaminarization 
is a great challenge. An objective analysis of the data yields no conclusive results indicating 
where and when relaminarization occurs. In fact after analyzing similar data sets for values 
of K between 1 x 10 -6 < K < 5.0 x 10“ 6 , it becomes apparent that relaminarization 
does not have a readily identifiable dependence on K. Most quantities change continuously 
with increasing K\ demonstrating no catastrophic changes when relaminarization is present. 
Therefore the need arises for a more suitable method of detecting relaminarization. 
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Figure 13: Boundary layer skin friction and friction velocity evolution for K = 4.10 x 10 6 
and K = 2.05 x 10~ 6 . 
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ce intensity evolution for K = 4.10x10 6 and K = 2.05x 10 -6 . 
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Figure 14: Maximum turbuleni 















3.2.3 Turbulent Burst-Sweep Event Detection 

As previously stated, relaminarization of a turbulent boundary layer is thought to be caused 
by the domination of the pressure forces over the slowly responding Reynolds stresses [3] ,[6], [10], 
but as Fernholz and Warnack[8] found, “this result excludes a parameter-connected criterion 
for the onset of relaminarization and for the breakdown of the standard logarithmic law\” 
These results are consistent with the current findings of this study. Indeed, this poses a 
serious problem if we are to assess the influence of relaminarization on any given engineering 
problem. 

In search for a better understanding of the flow physics associated with turbulent bound- 
ary layers exposed to large favorable pressure gradients, Fernholz and Warnack [7], [8] made 
detailed flow field measurements and documented the turbulent properties of the boundary 
layer. They found a gradual change in the turbulence properties of the turbulent boundary 
layer wdien exposed to large favorable pressure gradients; consistent with the claim that the 
relaminarization process is not a catastrophic event. In addition, Fernholz and Warnack[8] 
documented a significant change in the near-wall turbulence quantities for a relaminarizing 
turbulent boundary layer. Focusing on the fluctuating w T all shear stress, they found large 
changes in the low T er and higher moments, skewness, flatness, and the Reynolds normal stress 
in the near-w^all region (j/ + < 20). These significant changes in the near- wall region may be 
indicative of large changes in the turbulent burst-sw^eep cycle of events. If so, the turbulent 
bursting rate is a measurable quantity and may w T ell be a significant factor in determining 
the onset of relaminarization. 

Using Quadrant Splitting Analysis, a conditional sampling technique often used to detect 
burst-sw r eep events, the turbulent bursting frequency can be computed from time signals of 
fluctuating u and v velocity measured via an X-w'ire in the near-w T all region. The turbulent 
bursting frequency is defined as the time between ejections and can be nondimensionalized 
as 


n = ^ (is) 

using the freestream velocity and the boundary layer thickness. The higher the turbulent 
burst frequency, the larger the time between detected burst-sweep events, or, in other w r ords, 
the louver the turbulent burst rate. 

Applying Quadrant Splitting Analysis to the u and v time signals measured in the near- 
w r all region of the wind tunnel, the turbulent burst frequency is extracted and is plotted in 
Figure 16. For the Case 2 ( K = 2.05 x 10~ 6 ) boundary layer, there is a slight increase in 
the turbulent burst frequency through the contraction; this indicates a slight decrease in the 
turbulent burst rate as the flow r travels through the contraction and is most likely explained 
by the presence of the strong favorable pressure gradient acting as a stabilizing influence on 
the turbulence production. For the Case 1 (K — 4.10 x 10 -6 ) boundary layer, the turbulent 
burst frequency follows the Case 2 flow until approximately x = 3.94 m (155 in) w T here the 
burst frequency dramatically increases before relaxing to the exit values of Case 2. This 
increase in the burst frequency is associated with a large reduction in the turbulent bursting 
rate; a large reduction in turbulent bursting rate is consistent with relaminarization of the 
turbulent boundary layer. Tentatively, it seems the strong favorable pressure gradient acts 
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Figure 16: Turbulent burst frequency evolution calculated using Quadrant Splitting Analysis. 


to reduce the bursting rate in a turbulent boundary layer slightly, but imposed on this mild 
reduction in a turbulent boundary layer is an even larger reduction in the turbulent bursting 
rate when relaminarization is present. 

4 Relaminarization and the NASA Dryden FTF-II 

4.1 Design Concept 

The flight testing phase of the project is designed to determine the flow physics of turbulent 
boundary layers undergoing relaminarization at Reynolds and Mach numbers typical of high- 
lift systems found on commercial aircraft. In keeping with this focus, the flight test model on 
the FTF-II must approximate the pressure gradient and flow field found on the main element 
leading edge in these high-lift systems. The current design incorporates the upper half of 
the FTF-II with a horizontal splitter plate mated to the FTF-II. The horizontal splitter 
plate has a rounded leading edge which acts to accelerate the flow around the leading edge. 
To supplement this flow’ acceleration, a symmetric airfoil is mounted below’ the leading edge 
nose. The leading edge boundary layer is exposed to the favorable pressure gradient , w^hich 
is associated with the off-surface airfoil suction peak. It is proposed that the airfoil be 
mounted to a servo motor thereby allowing the airfoil to be pitched w^hile in-flight and 
increasing the number of tests that can be performed per flight. This design concept is 
illustrated in Figure 17. 

A preliminary analysis using a commercially available panel code has been performed on 
a representative model with a cylindrical leading edge for a small number of configurations. 
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Figure 17: FTF-II flight test schematic featuring horizontal splitter plate with rounded 
leading edge and rotatable symmetric airfoil for modifying pressure field. 



Figure 18: F-15B/FTF-II flight envelope with the current computational analysis location 
indicated. 
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All calculations were performed assuming a Mach number, Moo — 0-5, at an altitude of 
10,000 ft as indicated in Figure 18. This corresponds to a unit Reynolds number for the 
FTF-II of Re = 7.57 x IQ 6 . Two leading edge radii were chosen, R/L — 0.1 and R/L = 0.2, 
(where L = 2.69 m (106 in) is the length of the FTF-II). The chord of the symmetric airfoil 
was chosen to be c/L = 0.3 and the thickness to be t/L — 0.05 for both cases. The gap, 
defined as the distance between the leading edge of the main body (FTF-II) and the leading 
edge of the symmetric airfoil at zero angle of attack relative to the main body, was chosen 
to be gap/L = 0.2 and gap/L — 0.3 for the 10% and 20% radius case respectively. For 
both leading edge radii, computations were made at angles of attack of 0°, 5 3 , and 10° 
w r ith respect to the main body. Figure 19 shows a plot of the relaminarization parameter 
K vs. non-dimensional distance, z/L, for the R/L = 0.1 leading edge radius at the three 
angles of attack denoted previously. The value of K = 3 x 10“ 6 is indicated as well. 
Clearly Figure 19 indicates sufficiently large values of the relaminarization parameter to 
achieve relaminarization on the leading edge, if K is the appropriate governing parameter. 
Figure 20 is a similar plot for the R/L = 0.2 leading edge radius and it also indicates the 
ability to achieve relaminarization on the leading edge, if, again, K is indeed the appropriate 
parameter. Therefore, with a simple model configuration such as the one illustrated in 
Figure 17, sufficiently strong flow acceleration can be developed to achieve relaminarization 
on the FTF-II assuming K is the appropriate design variable. 
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K vs x/L 

FVL ■ 0.1, VI - 0.05, cA. « 0.3, gap/L * 0.2 



Figure 19: Computation #1 on FTF-II model leading edge indicating K values large enough 
to achieve relaminarization. 



Figure 20: Computation #2 on the FTF-II model leading edge indicating a K value large 
enough to achieve relaminarization. 
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Measured Re e at x = 3.4 m (132 in) 



Figure 21: Re# behavior for given acceleration parameter K. 

5 Recommendations 

5.1 Short Coming’s of the Relaminarization Test Facility 

While the relaminarization test facility has yielded valuable data, there are several concerns 
regarding the current approach. Chief among these is the unavoidable fact that to achieve 
increasing values of the acceleration parameter K , the initial Reynolds number must decrease, 
as can be seen in Figure 21. The problem exists because the larger K values are achieved 
at relatively small velocities, which makes sustaining a fully developed turbulent boundary 
layer difficult. In fact, there is evidence which suggests the presence of low Reynolds 
number effects in the initial turbulent boundary layer for the V ^ = 4 — and U^o = 5 ® 
cases. These effects can be seen in Figures 7-9 and 15. In addition, by decreasing the 
initial Reynolds number to achieve greater values in the acceleration parameter, the effect of 
pressure gradient becomes more and more dominant over a turbulent boundary layer whose 
turbulence is becoming weaker; making it difficult to make definitive statements about the 
effect of K on relaminarization. 

5.2 A New Relaminarization Facility 

Future work will focus on the design and construction of a new relaminarization facility in 
the University of Notre Dame’s Hessert Center’s 1.5m xl.5m (5ft x5ft) atmospheric wind 
tunnel. The reasons for moving from the current 0.6 m x 0.6 m (2 ft x 2 ft) indraft wind 
tunnel are threefold: 
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1.5 m x 1.5 m (5 ft x 5 ft) Atmospheric Wind Tunnel 



NOT TO SCALE 


Figure 22: Schematic drawing of proposed relaminarization test facility in 1.5 m x 1.5 m 
(5 ft x5ft) atmospheric wind tunnel. 

1. Operate with a fixed initial and higher Reg while maintaining a variable K environment 
through the use of a adjustable wall contour, 

2. Increase the initial Reynolds number to ensure a fully turbulent boundary layer, 

3. Increase the initial boundary layer thickness to allow for greater Tvall-normal spatial 
resolution in the measurements. 

In the current facility, to achieve greater values of the acceleration parameter, the 
initial Reynolds number must decrease. To address this problem, the new test facility will 
incorporate a adjustable well contour. Maintaining the linear contour used in the current 
0.6 m x0.6 m test facility to achieve a constant K environment , the 1.5 m x 1.5 m test facility's 
linear contour will be hinged at the ceiling and the angle of inclination to the top wind tunnel 
wall will be adjustable. Figure 22 is a schematic of the proposed relaminarization test facility 
in the 1.5m x 1.5 m atmospheric tunnel. Unlike the current test facility in which selection 
of the initial Reynolds number dictates the K value, adjustment of the wall contour in the 
1.5 m xl.5m test facility will provide the desired K value for any prescribed initial Reynolds 
number. 

It is planned to conduct the relaminarization experiments in the 1.5 m x 1.5 m atmospheric 
tunnel with an initial velocity of Uoo = 5 — , which corresponds to a initial Reynolds number 
of Re$ = 5660. Three constant values of the acceleration parameter will be studied: 
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1. mild favorable, ^ 1.5 x 10" 2 * * * 6 

2. moderate favorable, K ~ 3.5 x 10~ 6 

3. strong favorable, K ~ 5.0 x 10 -6 . 

In 1998, Fernholz and Warnack[8] published the most extensive work on turbulent bound- 
ary layers exposed to favorable pressure gradients in nearly twenty years. There work was 
conducted at initial Reynolds numbers of Re# = 862 and 2564 for K < 4.0 x 10“ 6 . A lit- 
erature search has found that nearly all previous studies were conducted at initial Reynolds 
numbers of Re^ 3000 A Thus this new facility will permit us to investigate the response of 

the turbulent boundary layer to favorable pressure gradients at much larger initial Reynolds 
numbers and greater flow acceleration than previously explored. 

Although larger initial Reynolds numbers are useful and needed, the true novelty of this 
experiment lies in the continuing use of a constant K region. Previous studies have relied 
upon a variable K region that peaked at a given value. This new facility will maintain a 
constant K region for a significant duration and makes for a clean experiment. 

5.3 FTF-II Flight Testing 

At present, there is insufficient tools available to immediately proceed with the continued 
development of a FTF-II flight experiment. The current experimental evidence suggests 
that K is not a reliable indicator of relaminarization. Further more, our ability to detect 
relaminarization has, at best, been shown in ground testing to be severely limited at present. 
Before a suitable FTF-II flight test design is completed and flown, several questions must be 
addressed, chief among these are: 

1. Under what conditions will relaminarization occur and what is the governing parame- 
ter? 

2. When relaminarization is present, what is the appropriate tool necessary to detect 

relaminarization? 

The planned wind tunnel testing in the Hessert Center’s 1.5 m x 1.5 m (5 ft x5ft) atmo- 
spheric wind tunnel is designed to answer these questions. Upcoming tests in the atmo- 
spheric wind tunnel will shed light on the nature of relaminarization and, it is believed, lead 
to the development of a new governing parameter. Once this is determined, the FTF-II 
design discussed previously may then be appropriately modified, so as to ensure relaminar- 
ization on the FTF-II in flight testing. 

Based on the experimental evidence to date, it is believed that turbulent burst-sw^eep 
measurements will be a viable detection criterion for determining w r hen and wriiere relami- 
narization occurs. Development of this diagnostic tool will continue in wind tunnel tests 
with the ultimate goal being to use flush mounted hot film arrays in flight testing and then 
applying a burst-sweep event detector algorithm to the data to determine, either in real-time 
or post-flight, the location and nature of relaminarization at Mach and Reynolds numbers 
more typical of commercial high-lift systems. 

^nly one case w*as found with a greater initial Reynolds number; Patel and Head [5] in 
1968 investigated flow" in a pipe with Re# = 5900. 


26 


6 Summary 

Turbulent boundary layer relaminarization has been shown to potentially exist on high- 
lift systems from the results of van Dam, et al.[2] and a possible mechanism by which 
relaminarization may manifest itself has been illustrated. Under the current research effort 
supported by the NASA Dry den Flight Research Center, ground testing facilities have been 
developed at the University of Notre Dame in order to further our understanding of turbulent 
boundary layer relaminarization and to assess the importance of relaminarization in high-lift 
systems. To date, testing has been conducted in the 0.6 m x 0.6 m (2 ft x 2 ft) subsonic indraft 
w T ind tunnel at the Hessert Center for a range of constant acceleration parameters. Results 
(shape factor evolution, skin friction evolution, maximum turbulent intensity evolution, log- 
law behavior, etc...) are consistent with the presence of relaminarization in the expected 
cases; although there is some question as to whether these are adequate discriminators for 
determining the presence of relaminarization. In addition, a preliminary concepts for an 
in-flight relaminarization experiment on the FTF-II has been explored. 

Experimental evidence suggests that the acceleration parameter, A", is not a viable pa- 
rameter for determining the onset of relaminarization. This poses a large problem, for 
designing a successful relaminarization flight test experiment, since there exists no defini- 
tive analysis tools to determine in advance when and wdiere relaminarization would occur. 
However, preliminary evidence supports the claim that the turbulent bursting rate may be 
a direct indicator of relaminarization. 

Future w T ork will focus on the development of a new relaminarization facility in the Uni- 
versity of Notre Dame’s Hessert Center’s atmospheric wind tunnel, w r hich will provide a clean 
experiment for the further study of relaminarization. It is intended to use this facility to fur- 
ther assess the viability of the acceleration parameter and to develop a new relaminarization 
parameter to aid in the development of the proposed FTF-II flight test experiment. In ad- 
dition, development w'ill continue on the use of the turbulent burst-sweep event as a method 
for detecting turbulent boundary layer relaminarization and applying a detector algorithm 
to flush mounted hot films in order to construct an in-flight relaminarization detection tool. 
This work will culminate in an in-flight relaminarization experiment for the FTF-II and the 
suitable flow diagnostic tools to detect relaminarization in flight. 
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